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Abstract 

Rift Valley fever is a vector-borne disease which causes high morbidity and mortality in live- 
stock and also can infect humans. In the event Rift Valley fever virus is introduced to the United 
States, mitigation strategies need to be rapidly set in place to prevent the disease from becom- 
ing endemic. Presented here is a stochastic network-based mathematical model of Rift Valley 
fever. Given a lack of empirical data for this region, this difference equation model uses stochas- 
tic parameters following several Pert distributions. A large directed, asymmetric network was 
established with 3621 nodes in the Uvalde and Panhandle areas of Texas, U.S.A. noting that the 
climate in Texas can be similar to that of the Rift Valley in Eastern Africa. This model predicts 
how Rift Valley fever can spread spatially based on biotic and antibiotic factors and identifies the 
relevant parameters in RVFV propagation. Humans, cattle, Aedes mosquitoes, Culex mosquitoes 
are independently treated with diff'erent contact networks to assess virus propagation and spread 
under diff'erent scenarios. The importance of the factors aff'ecting Rift Valley fever virus spread 
was assessed under diff'erent scenarios. Spatial eff'ects and climate factors are also addressed in 
the model. It has been shown that cattle movement, mosquito population and favorable climate 
for mosquito development contribute to the spatial spread of Rift Valley fever virus. 
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1. Introduction 

Rift Valley fever (RVF) was first identified in Egypt in 1931 [6] and is endemic in east and 
South Africa [9]. It is a mosquito-borne viral disease which mainly aff'ects ruminants and humans 
[17]. Rift Valley fever causes abortion in adult and death in newborn livestock [17]. Sheep, goats 
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and cattle are the most important affected animal hosts from an economic standpoint [9]. Humans 
and domestic livestock are the main hosts [5]. The virus is maintained by vertical transmission 
within Aedes mosquitoes [16]. High levels RVF activity is usually related to persistent, above 
normal rainfall and ENSO event in East Africa [15]. 

Ruminants usually acquire the infection from the bites of mosquitoes [18]. Aedes and Culex 
mosquitoes are thought of as the main vectors with respect to vector competence [5]. Aedes 
mosquitoes lay eggs in dry mud [30] and the eggs can survive for a long period [9]. After 
flooding, the eggs infected by Rift Valley fever virus (RVFV) can develop into infected adult 
mosquitoes [9]. 

Movement of infected animals is thought of as a possible pathway for introducing RVFV to 
a disease-free area [19]. In Egypt in 1977 [25], the trade of sheep from east Africa to Mecca in 
Saudi Arabia or in Yemen during a festival period was considered to be a likely pathway for the 
introduction of RVFV [2, 7]. As a consequence, a ban on livestock trade can be an efl'ective way 
to prevent further spread of RVFV [7] . 

Humans can acquire the infection from the bites of mosquitoes or bodily fluids of animals 
during butchering or slaughtering [1]. However, most human infections result from contact with 
animals [9]. Individuals working with animals, such as farmers and veterinarians, are the most 
vulnerable to RVFV infection [23]. Rift Valley fever virus infection causes severe influenza-like 
disease in humans with serious consequences such as blindness, or even death [17]. It has been 
reported that more than 200 humans died of RVF in Mauritania in 1987 [13]. There were 738 
reported human cases in Sudan including 230 deaths, in 2007 - 2008 [28]. It is possible that the 
number of human cases has been underreported, especially in rural areas [5]. Rift Valley fever 
virus has already spread to Yemen and Saudi Arabia [5] and the French island of Mayotte with 
multiple human cases reported [26]. Rift Valley fever virus may be introduced to the United 
States, similar to West Nile virus [10], and the consequences of RVFV transmission to humans 
in new locations should be focused on. 

Rift Valley fever models have been proposed; however, no model addresses spatial trans- 
mission during an outbreak [19]. Very few mathematical dynamic transmission models have 
explored mechanisms of RVFV circulation [19], although one dealt with the varied prevalence 
of RVF and mosquito abundance data from 1991 to 1996 [4]. A theoretical model in a closed 
system including Aedes and Culex mosquitoes and livestock population has been proposed by 
Gafl' et al. [11]. The authors suggest that RVFV can exist if the contact rate between vectors and 
hosts is high based on sensitivity analysis of the theoretical model parameters [19]. However, 
neither human hosts nor spatial efl'ects are considered in their model. A theoretical mathematical 
model on RVF with human hosts has been proposed [20] . The important parameters to the trans- 
mission of RVFV have been assessed in their model using sensitivity analysis. Human hosts are 
considered in their model, but spatial efl'ects and climate factors are not considered. For the Ferlo 
area of Senegal, a pond-level meta-population model was built assuming that Aedes mosquitoes 
are the only vector and rainfall is the only driving force [9] . It has been shown that within-Ferlo 
virus is possible to persist if the livestock move between ponds and the rainfall does not occur 
in all ponds simultaneously [9]. Although various transmission scenarios have been explored, 
no hosts were considered in the model. No hosts are considered in the model. A determinis- 
tic meta-population model including human hosts has been proposed by our group [29] and a 
study was carried out for South Africa. The model is based on RVFV spatial transmission during 
an outbreak, where a small network with three nodes has been established [29] . In reality, the 
propagation of the epidemic is very complicated, so a deterministic model is not always enough 
to reflect the dynamic spread of the disease. The real network can be large, complicated and 
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hierarchical. 

Proposed here is a stochastic network-based mathematical model with four competent species, 
that is, Aedes mosquitoes, Culex mosquitoes, cattle and humans. The incidence of infection can 
be predicted and mitigation strategies can be assessed with this model. Texas is roughly climati- 
cally similar to areas of the Rift Valley in eastern Africa where RVF is endemic . It's hypothesized 
that RVFV can be introduced to Texas based on empirical animal and movement data available 
for two regions in Texas. The nodes were chosen as 3526 farms in nine counties in South Texas: 
Bandera, Edwards, Frio, Kinney, Maverick, Medina, Real, Uvalde, and Zavala Counties [27], 
three livestock markets located in three counties: Uvalde, Medina and Frio counties [27], and 
92 feedlots located in the northern Panhandle area [27]. In all, 3621 nodes populated the large 
network. Different scenarios were then simulated to analyze how RVFV will spread and which 
factors that are playing an important role in RVFV spreading if infected mosquitoes arrive in the 
United States. The results have shown that cattle movement, mosquito population, and climate 
play an important role in the spatial spread of RVFV. 

2. The stochastic mathematical model 

Three different sets of nodes were considered in the model, namely, 3526 cattle farms, or 
ranches and 3 livestock markets in the southwest of Texas and 92 feedlots in the northwest part of 
the state [27]. The nodes were considered homogeneous within each set. This meta-population 
model is used to study RVF dynamics in the 3621 locations. Each farm, market or feedlot is 
represented as location /, /=!,•••, 3621. All the transitions to be discussed below are for 
location / at day t. 

2.1. The Difference equation model 

The model is based on the following assumptions. 

* Each farm, livestock market, feedlot is a node. The network is directed and asymmetric 
considering movement in different directions and rates of moving in and out of a node. 
Farms, feedlots and markets are in different categories respectively. 

* The animals include only cattle. 

* Once an animal leaves a node, it never returns to that node. 

* All animals are homogeneously mixed within a node. 

* Only RVFV transmission in farms and in feedlots is considered because markets are hosts- 
transit nodes. 

* The movement rate of cattle in infected compartment is assumed to be half of that of the 
cattle in other compartment in the absence of further information. Any realistic movement 
rate can be adopted in the model. 

* The mosquitoes will only diffuse to adjacent nodes in a time step. The destination of cattle 
movement is randomly selected and independent of distance. 

* The birth rate and death rate is in markets, birth rate is in feedlots. 

* Human birth, death, mortality and mobility are ignored. 
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Birth, death, infection, transfer between compartments, mosquitoes diffusion, and cattle 
movement occur at the end of each time step. 



2.1.1. Aedes mosquito vectors 

Below is the set of diffence equations for Aedes mosquitoes. The parameters and variables 
are described in 2.2. 

Pu{t + 1) - Pu{t) = {bi [Nu{t) - qihiit)] - OiPuit)) * 6A{t) (1) 
Quit + 1) - Quit) = (biqihiit) - OiQuit) * dA{t) (2) 

n n 

Su{t+\)-Su{t) = e,Pu{t)+ ^)iSijit)- Yj ^]jSu{t)-d,Su{t) 

j=ij*i j=Uj*i 
-p2lSuif)hiit)IN2i{t) (3) 

n n 

Eu(t + I) - Eu(t) = Yj ^-^lyW- Tj <^hEu(t) - diEu(t) 

+ /32i(t)S xiit)hi{f)IN2iit) - eiEuit) (4) 

n n 

hiit+\)-hiit)= Y 4 Y <^lMt) + OiQuit) - dihiit) + eiEuit) (5) 

Nu(t +l) = S u(t + 1) + Eu(t + 1) + hiit + 1) (6) 
Qu(t+l) = Qu(TQ[t/TQ\ + l) (7) 

A modified S-E-I model includes compartments P and Q which are uninfected eggs and in- 
fected eggs respectively [29]. There are biNuit) eggs laid, including biqiluit) infected eggs and 
biNiiit) -biqihiit) uninfected eggs each day. After the development period, OiPuit) uninfected 
eggs develop into susceptible Sidult Aedes mosquitoes and 6i Quit) infected eggs develop into in- 
fected Sidult Aedes mosquitoes every Tq days. The percentage of Aedes mosquitoes moving from 
location / to location j is 61;!^.. The change in the number of Aedes mosquitoes due to mobility in 
compartment X is given as Yj%ij^i ^j/^ijCO - ^Ij^uO)- The number of infected Aedes 

mosquitoes from contact with livestock is yS2i^ 1/(0^2/(0/^2/(0- Following the incubation period, 
£iEii(t) Aedes mosquitoes transfer from exposed state to infected state. Aedes mosquitoes dying 
naturally in compartment X is given as diXuit). Dead mosquitoes are removed from the com- 
partments. If the average temperature is above 60 and the average rainfall is above mm 
during the Tq days, then ^^(0 = 1. otherwise SaO) = 0. If the average temperature or rainfall 
during the period does not reach the threshold, new eggs are not produced and the infected eggs 
will not develop into adult mosquitoes. 

2.1.2. Culex mosquito vectors 



P^i(t + 1) - ^3/(0 = b^(Temp, precipitation, T^, t)N^i(t) - 6^(1 emp, t)P^i(t) (8) 

n n 

Sii(t+l)-Sii(t) = 0i(temp,t)Pii(t)+ Y ^^3/0- Y ^lSii(t)-d3Sii(t) 

-/323S3imi(t)/N2i(t) (9) 
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Eii(t + I) - Eiiit) = Yj ^l^ijit)- Yj o?ijE3i{t)-SiEiiit)-diEiiit) 

j=Uj*i j=lj*i 

+ /323S3imi(t)/N2i(t) (10) 

n n 

hi(t^l)-hi(t)= Yj ^]Mt)- Yj Ojlhi(t)^S3E3i(t)-d3l3i(t) (11) 

N3t(t + 1) = ^3/^ + 1) + E3t(t + 1) + ht(t + 1) (12) 

Above is a modified S-E-I model for Culex mosquitoes. Only uninfected eggs are included since 
the female Culex mosquito has not been demonstrated to transmit RVFV to its own eggs [11]. 
There are b^N^iit) eggs laid each day. After the development period, O^P^iit) eggs develop into 
susceptible adult Culex mosquitoes. The percentage of Culex mosquitoes moving from location 
/ to location j is oj^j. The change in the number of Culex mosquitoes due to movement in com- 
partment X is given as j^i <^\^?>j(f) - Z}=i j^i ^^^^3/(0- The number of Culex mosquitoes 
acquiring infection from livestock is denoted by I323S 3i(t)l2i(t)/N2i(t). After the incubation pe- 
riod, s^E^iit) Culex mosquitoes transfer to the infected compartment. The Culex mosquitoes 
removed from the compartment X due to natural death is d^X^iit). The total number of Culex 
mosquitoes is also variable and equals the sum of the Culex mosquitoes in all the compartments. 

2.1.3. Livestock hosts 



S2i{t+\)-S2i{t) = b2{t)d,ii)N2i{t)+ Yj 4^2/0- Y ^p2i(t)-d2SAi)S2im2i(t)/K2 

-Pi2S2i{t)Iu(t)INu{t)-l332S2imi{t)IN3iit) (13) 

n n 

E2i(t + I) - E2i(t) = 4^2/0- J] <olE2i(t)-d2dd(i)E2i(t)N2i(t)/K2-S2E2i(t) 

+ PnS2i{t)hi{t)INu{t) +/332S2imi(t)/N,i(t) (14) 

n n 

hiit + I) - hiiO = Yj "^^jMt)- Y (^lhi{t)-d26d{i)hi{t)N2i{t)IK2 + e2E2i{t) 

-y2hiit)-^l2hi(t) (15) 

n n 

R2iit+\)-R2iit)= Y 4^2,(0- Y ^l^2i(t)+72hi(t)-d2Sdm2im2i(t)/K2 (16) 

N2i(t + 1) = S2i(t + 1) + E2i(t + 1) + hiit + 1) + R2i(t + 1) (17) 

It's assumed that cattle can be transported from farms to markets or feedlots, or from markets 
to farms or feedlots, but not from feedlots to farms. The total number of newborn livestock in 
location / is (0^2/(0 each day. The percentage of cattle moving from location / to location 
j is ojfj. The change in the number of livestock in compartment X due to mobility is given as 

I.%ij^i^jiX^^\t) - I,%ij^iOjlX^^\t). After the incubation period, ^2^2/(0 livestock transfer 
from exposed state to infected state. The number of cattle infected by Aedes mosquitoes and 
Culex mosquitoes is denoted as jSuS 2i(t)hi(t)/Nii(t) and yS325' 2/(0^3; (0/^3/(0. respectively. After 
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the infectious period, jihiit) cattle recover from RVFV infection. The number of dead cattle in 
each compartment X is given as ^2^2/^^2/(0/^2 in which K2 is the carrying capacity of cattle in 
each node. There are jdihiit) cattle dying of RVFV infection. The variables di){i) and dd{i) are 
used to differentiate different type of nodes. If location / is a farm, then Si^ii) = 1, Sd(i) = 1. If 
location / is a market, then Sh(i) = 0, ^^(0 = 0. If location / is a feedlot, then ^^(0 = 0. 

2.7.4. Human hosts 



S4i(t + 1) - ^4/(0 = -PuSAi{t)hi{t)INu{t) -fi24S4i(t)l2i(t)/N2i(t) - I^^aS Ai{t)hi{t) IN^i{t) (18) 

E,iit + 1) - E.iit) =/3uS4imi(t)/Nu(t) +P2,fiSM{t)hi{t)IN2i{t) + PmS ,i{t)hi{t)IN^i{t) 



Assume that the susceptible human population only includes those working with animals, such 
as veterinarians and farmers, because they are at higher risk of RVFV exposure and infection 
[23]. The number of humans infected by Aedes mosquitoes, Culex mosquitoes, and livestock 
is ySi4^4/(0/i/(0/A^i/(0, P24S4i(t)l2j(t)/N2i(t), SiYid P34S 4i(t)hi(t) /Nsiit) respectively. There are 
£:4£'4/(0 humans transferring to the infected compartment after the incubation period and 74/4/(0 
humans recover from RVFV infection. 

2.2. Descriptions of parameters and variables in the model 

The descriptions of the parameters of the following equations are listed in Table 1 . 
The variables in the equations are as following: 
Pi/(0=The number of uninfected A ^J^^ mosquito eggs in location / 
2i/(0=The number of infected A ^J^^ mosquito eggs in location / 
S i/(0=The number of susceptible Aedes mosquitoes in location / 
£'i/(0=The number of exposed A ^J^^ mosquitoes in location / 
/i/(0=The number of infected A ^J^^ mosquitoes in location / 
A/^i/(0=The total number of Aedes mosquitoes in location / 
^2/(0=The number of susceptible livestock in location / 
£'2/(0=The number of exposed livestock in location / 
/2/(0=The number of infected livestock in location / 
^2/(0=The number of recovered livestock in location / 
A^2i (0=The total number of livestock in location / 
P3/(0=The number of uninfected Culex mosquito eggs in location / 
5'3/(0=The number of susceptible Culex mosquitoes in location / 
£'3/(0=The number of exposed Culex mosquitoes in location / 
/3/(0=The number of infected Culex mosquitoes in location / 
A^3/(0=The total number of Culex in location / 
5'4/(0=The number of susceptible humans in location / 
£'4/(0=The number of exposed humans in location / 
/4/(0=The number of infected humans in location / 
7?4/(0=The number of recovered humans in location / 



Uiit + 1) - /4/(0 
R4i{t+\)-R4i{t) 



S4E4i(t) 

S4E4i(t) - 74^4/(0 
74^4/(0 



(19) 
(20) 
(21) 
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3. The network 



The movement between the three sets of nodes is shown in Figure 2(a). The network is 
directed and asymmetric. The distance between any two farms is compared with an assumed 
radius which is defined as 2 km here. If the distance is smaller than the radius, then the two farms 
are neighbors for mosquito movement and there is a link between them. An example of the 
mosquito movement network is shown in Figure 2(b). All nodes are shown in Figure 3(a), where 
links are omitted due to their large number. An example of the links of a single farm is shown in 
Figure 3(b). The exact location of each farm, market and feedlot is obscured because those data 
are not available to public. The nodes are uniformly distributed in each county according to the 
total number of farms in each county [27]. Some markets are randomly assigned as neighbors of 
each farm and some farms are randomly assigned as neighbors of each market. The neighbors 
are the destination of the cattle. Farms and markets are thought of as sources of cattle to feedlots. 
The destination, markets or feedlots, from cattle farm is independent of distance as are cattle 
from markets. The movement rate is inverse to distance. The nodes are distributed in each 
county according to the total number of feedlots in each county in the Panhandle area [27]. The 
initial number of susceptible cattle in each farm or feedlot for numerical simulation is assumed 
according to the total number of cattle in each county [27] . The color of the node is determined 
by the number of infected cattle in each farm, market or feedlot. Green represents the susceptible 
cattle and humans respectively. Initially, there are no infected nodes in any of the counties. 

4. Case Study: Texas, U.S.A. 2011 

Case study was carried out for Texas, U.S.A from January, 2005 to December, 2006. The 
parameters of the above equations are listed in Table 1 . The parameters with given ranges are 
variables following severals Pert distributions. Uniform distribution has equal probability for any 
value within the range. Pert distribution has a smooth trend and flat tail as Gaussian distribution 
as well as a flexible most possible value. Therefore, Pert distribution is chosen for our preference. 
The movement rates of livestock are shown in Table 2. The movement rate is assumed to be the 
maximum percentage of animals moving [8]. The number of farms and cattle in each of the nine 
counties [21] is shown in Table 3. 

where: 

^m(0=the number of markets connected to farm / 
^/(/)=the number of farms connected to market / 
^//e(/)=the number of feedlots connected to farm / 
^m/e(0=the number of feedlots connected to market / 
diSij=\hQ distance between node / and node j 

The difl'usion rates of Aedes and Culex mosquitoes are shown below [24]. 

I 3/0 if the nodes are disjoint 

^ij ~ ^ij ~ I diff/dfj if two nodes share a border 

where dijis the distance between the centers of node / and node j [24] and diff is a difl'usion like 
parameter within the range (830, 8300)m^/day [24]. 

The development rate and egg laying rate of Culex mosquitoes [12] are shown in Equation 
(22) and Equation (23). The parameters adopted are derived from data concerning West Nile 
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Figure 1 : The relationship between Z73 , 63 and climate 
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in Uvalde study area 



(a) Animal movement location (b) An example of mosquito network 

Figure 2: Movement diagrams 




virus in 2010 in the northern US [24]. This is the best model currently available for us at the mo- 
ment. Any more realistic models or parameters can later be adopted in the model once available. 
The egg laying rate changes with moisture [12] and development rate changes with temperature 
[12] are shown in Figure 1(a) and Figure 1(b) , respectively. The egg laying rate and development 
rate of Culex mosquitoes in Texas from January, 2005 to December, 2006 are shown in Figure 
1(c) and Figure 1(d), respectively. The development rate in Figure 1(d) is enlarged 100 times 
to be seen clearly. As is shown, the development rate of the Culex mosquitoes changes with 
temperature. The colder it is, the fewer Culex eggs are produced. If the temperature is too low, 
the eggs won't develop into adult mosquitoes. If the temperature is too high, the lifespan of the 
mosquitoes is shortened so that the development rate decreases. Moisture index is the difference 
between precipitation and evaporation as shown in Equation (24). The lower the moisture index, 
the fewer adult mosquitoes are produced because low moisture index represents low precipita- 
tion and high evaporation. For some days, the missing precipitation data from January, 2005 to 
December, 2006 [22] missing are assumed to be zero. The evaporation data are calculated using 
the formula in Equation (25) [14]. One set of climate data was applied to the Uvalde study area 
and another set of climate data was applied to Panhandle area. The rainfall and temperature de- 
termine the presence of the mosquitoes. They are necessary for mosquitoes to survive and breed. 
The presence of the infected eggs is clearly important to begin the simulated epidemics. 



298.15 i+,^^[^*(^_^_L_)] 

Emax 

b^iTemp, precipitation, T,, t) = bo ^ — ^ Moisture(t)-Emean - (23) 

i + ex PI ^^^^ J 

t 

Moistureit) = ^ precipitation(D) - Evaporation{D) (24) 
100(Temp(t) -h 0.006/z)/(100 - latitude) 



Evaporation{t) 



80 - Tempit) 
l5(Temp(t) - TM) 



80 - Tempit) 



mm I day (25) 



Where: 

Tempit) =air temperature in units of ^C [14] 
latitude =the latitude (degrees) [14] 
Tdit) =the mean dew-point in units of [14] 
h =the elevation (meters) [14] 

Various scenarios have been examined and listed in Table 4. Simulations were run with 200 
different sets of parameters using a daily time step. The climate data from January, 2005 to 
December, 2006 [22] in Texas was adapted into the model. The spread of RVFV with time is 
shown in Figure 4. The parameter Tq is set to 15 as an example. Digital numbers are used 
to distinguish farms. As an illustration, the distance between Farm 3526 and other 4 farms is 
less than 2 km. So Farm 3526 has 4 neighbors. The change in the number of infected cattle or 
humans in farms during an outbreak is the main focus of this case study. Therefore, the number 
of infected cattle and humans in feedlots are not illustrated in the following figures. Assume 10 

10 



infected Aedes mosquito eggs are introduced to Farm 3526. The number of infected humans and 
cattle in Farm 3526 and its neighbor farms is shown in Figure 4(a) and Figure 4(b). To see the 
number of infected cattle and humans more clearly, the infected humans and animals only in the 
neighbor farms of Farm 3526 is shown in Figure 4(c) and Figure 4(d). If there are no infected 
eggs in all nodes at the beginning, then the percentage of infected animals and humans is zero in 
all farms. Farm 3440 is one of the neighbors of Farm 3526. The percentage of infected cattle and 
humans in the neighbor nodes of Farm 3440 shown in Figure 4(e) and Figure 4(f). The infected 
cattle and humans also arise in the neighbors of Farm 3440. If there is no movement, then there 
are no infected cattle or humans in the neighbors of Farm 3526 as shown in Figure 4(g) and Figure 
4(h). The spatial spread of RVFV with 10000 and 10 eggs in a small farm of Bandera county 
initially is shown in Figure 5. The percentage of the infected cattle and humans reaches the peak 
in March for Farm 3526 and around May for its neighbors. Rift Valley fever virus spreads much 
more quickly with more infected eggs initially as shown in the figure. If 10 infected eggs only 
exist in a farm with very few cattle or humans at the beginning, then cases likely won't spread 
quickly to other counties as shown in Figure 5(a) and Figure 5(b). Therefore, the infection of a 
small farm with few infected eggs does not aff'ect other nodes like an epidemic in a major node. 
The spatial spread of RVFV with 10 infected eggs in a large farm of Zavala county initially is 
shown in Figure 6. The percentage of infected cattle and humans on the lOOth day is shown in 
Figure 6(a) and Figure 6(b) respectively. By way of more comparison of Figure 5 and Figure 
6, Rift Valley fever virus spreads more quickly with infected eggs in a larger farm initially. The 
percentage of infected cattle and humans on the lOOth day with infected eggs in 9 large farms of 
9 counties at the beginning is shown in Figure 7(a) and Figure 7(b). RVFV has spread to other 
counties much more quickly as shown in the figures because more counties have infected eggs at 
the beginning. The infected cattle and humans also arise in the nodes of other counties because of 
livestock movement. If infected eggs are in large farms of diff'erent counties, RVFV will spread 
quickly. If 10 infected eggs are only in Farm 3526 at the beginning and there is only mosquito 
movement, the percentage of infected cattle and humans on the lOOth day is shown in Figure 8. 
If 10 infected eggs are only in Farm 3526 at the beginning and there is only cattle movement, 
the percentage of infected cattle and humans on the lOOth day is shown in Figure 9. The spatial 
spread of RVFV is more rapid with cattle movement than with only mosquito movement. The 
percentage of infected cattle and humans being infected with mosquito movement only, and 10 
and 10000 infected eggs initially, are shown in Figures 10(a), 10(b), 10(c), 10(d) and Figure 
10(e), 10(f), 10(g), 10(h), respectively. As is shown in Figure 10, the number of infected eggs 
greatly aff'ects the percentage of infected cattle and humans. The percentage of infected cattle and 
humans changing with time with only cattle movement and with diff'erent numbers of infected 
eggs are shown in Figure 1 1 . Therefore, the movement of cattle contributes to the spatial spread 
of RVFV with the given initial condition and movement rate as shown in Figure 8, Figure 9, 
Figure 10 and Figure 11. This gives us some information on potential for mitigation of RVFV 
spread. Banning animal trade and movement may help prevent further spread of RVFV during 
an outbreak. 

5. Discussions 

Presented here is a mathematical model with stochastic parameters describing the spatial 
spread of RVFV during an outbreak. The network is a large hierarchical network with three dif- 
ferent sets of nodes which are farms, markets and feedlots. The nodes are homogeneous within 



11 



Parameter 


Description 


Range 


Most 

possible 

value 


Units 


Source 


iZ 


adequate contact rate: Aedes to livestock 


(0.0021,0.2762) 


0.1392 


1/day 


[11] 




adequate contact rate: livestock to Aedes 


(0.0021,0.2429) 


0.1225 


1/day 


[11] 




adequate contact rate: livestock to Culex 


(0.0000,0.3200) 


0.16 


1/day 


[11] 


13^2 


adequate contact rate: Culex to livestock 


(0.0000, 0.096) 


0.04 


1/day 


[11] 


Pu 


adequate contact rate: Aedes to humans 


(0.001,0.002) 


0.0015 


1/day 


Assume 




adequate contact rate: livestock to humans 


(0.00004, 0.00008) 


0.0006 


1/day 


Assume 




adequate contact rate: Culex to humans 


(0.0005,0.001) 


0.00525 


1/day 


Assume 




recover rate in livestock 


(2,5) 


3.5 


1/day 


[11] 


1/74 


recover period in humans 


(4,7) 


5.5 


1/day 


[3] 


1/di 


longevity of Aedes mosquitoes 


(3,60) 


31.5 


days 


[11] 


l/d2 


longevity of livestock 


(360, 3600) 


1980 


days 


[11] 


I/J3 


longevity of Culex mosquitoes 


(3,60) 


31.5 


days 


[11] 


bi 


number of Aedes eggs laid per day 


1/(3,60) 


1/31.5 


1/day 


[11] 


b2 


daily birthrate of livestock 


1/(360, 3600) 


1/1980 


1/day 


[11] 


Z?3 


number of Culex eggs laid per day 


1/(3,60) 


1/31.5 


1/day 


[11] 


1/^1 


incubation period in Aedes mosquitoes 


(4, 8) 


6 


days 


[11] 


1/62 


incubation period in livestock 


(2,6) 


4 


days 


[11] 


1/63 


incubation period in Culex mosquitoes 


(4, 8) 


6 


days 


[11] 


1/64 


incubation period in humans 


(2,6) 


4 


days 


[20] 


U2 


mortality rate in sheep 


(0.025,0.1) 


0.0375 


1/day 


[11] 


Cl\ 


transovarial transmission rate in livestock 


(0,0.1) 


0.05 


1/day 


[11] 


IM 


development time of Aedes 


(5, 15) 


10 


days 


[11] 


1/^3 


development time of Culex 






days 


[11] 


^2 


carrying capacity of livestock 


100000 






Assume 




total number of farms 


3526 






[21] 




total number of feedlots 


92 






[21] 




total number of markets 


3 






[21] 


Tq 


development cycle of infected Aedes eggs 










A 


parameter of the development rate 


0.25 






[12] 


HA 


parameter of the development rate 


28094 






[12] 


HH 


parameter of the development rate 


35692 






[12] 


TH 


parameter of the development rate 


298.6 






[12] 


bo 


minimum constant fecundity rate 





1/day 




[12] 


Emax 


maximum daily egg laying rate 


20 


1/day 




[12] 


Emean 


the mean of the daily egg laying rate 





1/day 




[12] 


Evar 


variance of function 


12 


1/day 




[12] 


K 


Kelvin parameter 


273.15 


Kelvin 




[12] 



Table 1 : Parameter ranges for numerical simulation 
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i 


j 


range 


source 


farm 


market 


6Qn%linm{i)^dij) 


[8] 


market 


farm 


60.7% l{nf{i) ^ dij) 


[8] 


farm 


feedlot 


10.9% /(nffe(i)^dij) 


[8] 


market 


feedlot 


10.9% /(n^feii) ^ dij) 


[8] 


feedlot 


farm 





[8] 


feedlot 


market 





[8] 



Table 2: Movement rate of livestock oj^j 



County 


Cattle and calves 


cattle farms 


source 


Bandera 


11,254 


496 


[27] 


Edwards 


22, 247 


264 


[27] 


Frio 


51,411 


471 


[27] 


Kinney 


11,672 


125 


[27] 


Maverick 


25,135 


183 


[27] 


Medina 


55,759 


1364 


[27] 


Real 


4,830 


133 


[27] 


Uvalde 


52, 366 


316 


[27] 


Zavala 


66,641 


174 


[27] 


Total 




3526 


[27] 



Table 3: The number of farms and cattle in each county of Uvalde study area in 2007 



Scenarios 


Effect 


No infected Aedes eggs in any nodes 


No infected cattle or humans in any nodes 


Only 10 infected eggs are in a large farm in 
Zavala county at the beginning 


Infected cattle or humans appeared in other 
counties 


Only 10 infected eggs are in farm 1 (a small 
farm in Bandera county) at the beginning 


No infected cattle or humans observed in its 
neighbors. 


The infected eggs are in 9 farms of 9 differ- 
ent counties at the beginning 


Infected cattle or humans appeared in other 
counties 


Initially infected eggs are in a large farm 
in Zavala county and considering mosquito 
diffusion only 


Infected cattle and humans spread to other 
nodes 


Initially infected eggs are in a large farm in 
Zavala county and considering cattle move- 
ment only 


Infected cattle and humans spread to other 
nodes 


No cattle movement or mosquito diffusion 
and infected eggs are in Farm 3526 at the 
beginning 


No spatial spread of infected cattle and hu- 
mans 



Table 4: Scenarios with the climate selected during the most favorable period for mosquito development 
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100 200 300 400 

Time(Day) 

(a) Percentage of infected cattle in Farm 
3526 and its neighbors 




100 200 300 400 

Time(Day) 

(b) Percentage of infected humans in 
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100 200 300 400 

Time(Day) 

(c) Percentage of infected cattle in the 
neighbors of Farm 3526 
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100 200 300 400 

Time(Day) 



(g) Percentage of infected cattle in Farm 
3526 and its neighbors 




Tinne(Day) 

(h) Percentage of infected humans in 
Farm 3526 and its neighbors 



Figure 4: Numerical simulation results with Tq = 15 where the values are averaged over 200 realizations 



each set. The spread of RVFV is estimated not only within farms but also between farms, mar- 
kets, and feedlots. Evaluating several likely scenarios, the most important factors are identified. 
Model results indicate 

* Local transmission and trans-location transmission are important in accelerating the spread 
of RVFV. 

* The mosquito population and favorable climate for mosquito development contributes to 
the spread of RVFV among farms. 

* The route of moving from farms to livestock markets and subsequently to diff'erent farms 
expands the area of infection. 

The existence of infected eggs in large farms contribute to the rapid spread of RVFV. The model 
is a general model into which parameters in any location can be easily adapted. Climate factors 
are also involved and analyzed. The diff'usion of mosquitoes and the trans-location movement of 
livestock contribute to the spatial spread of RVFV. Starting with infected eggs in one large node, 
the infection of RVFV arises in more and more nodes later on. The infected Aedes mosquitoes 
diff'use to adjacent farms and hatch infected off'spring. The cattle can be infected by being fed 
on by Aedes mosquitoes. Culex mosquitoes amplify the infection following blood meals on 
infected cattle. The infected cattle transported to markets or feedlots are seeds of the infection in 
markets or feedlots. Similarly, the virus can be introduced to farms by purchasing infected cattle 
from markets. The parameters for the development rate and egg laying rate are not uniquely 
measured for Texas during 2005 to 2006. Rift Valley fever virus can spread quickly to the whole 
network once it is introduced into a large node. The actual animal movement rate from farms 
to markets, from markets to farms, from farms to feedlots, from markets to feedlots are not 
available to us due to confidentiality. Therefore, the incidence data can be inaccurate; however, 
the model gives us an idea of the pathways for rapid spread of RVFV among farms and counties. 
It shows that livestock trade and movements play an important role in RVFV spread between 
counties, restricting livestock trade and the movement of cattle during a RVF outbreak will 
reduce transmission. 
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0-0.000001 
0.000002-0.00001 
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(a) Percentage of infected cattle on the 100th 
day 
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(c) Percentage of infected cattle on the 100th (d) Percentage of infected humans on the 
day 100th day 



Figure 5: The dynamics of nodes in Uvalde study area with 10 and 10000 infected eggs in a small farm of Bandera 
county respectively at the beginning of the simulation 
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(a) Percentage of infected cattle on the lOOth 
day 
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Figure 6: The dynamics of nodes in Uvalde area with 10 infected eggs in a single large farm of Zavala county initially 
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100th day 



Figure 7: The dynamics of nodes in Uvalde area with 10 infected eggs in 9 large farms of 9 counties initially 
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0-0.000001 
0.000002-0.000010 
■ 0.000011-0.627369 

(a) Percentage of infected cattle on the lOOth 
day 



■ 0-0.000001 

0.000002-0.000010 
- 0.000011-0.032727 

(b) Percentage of infected humans the lOOth 
day 



Figure 8: The dynamics of nodes in Uvalde area with only mosquito diffusion and 10 infected eggs in a large farm of 
Zavala county initially 




0-0.000132 ■ 0-0.000002 

0.000133-0.002365 0.000003-0.000020 
■ 0.002366-0.435588 ■ 0.000021-0.035700 

(a) Percentage of infected cattle on the lOOth (b) Percentage of infected humans on the 

day lOOth day 

Figure 9: The dynamics of nodes in Uvalde area with only cattle movement and infected eggs in a large farm of Zavala 
county initially 
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neighbors of Farm 3526 with 10 infected 
eggs in Farm 3526 of Zavala county at 
the beginning considering mosquito dif- 
fusion only 
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(e) Percentage of infected cattle in Farm 
3526 and its neighbors with 10000 in- 
fected eggs in a large farm of Zavala 
county at the beginning considering 
mosquito diffusion only 




Time(Day) 



(f) Percentage of infected humans in 
Farm 3526 and its neighbors with 10000 
infected eggs in a large farm of Zavala 
county at the beginning considering 
mosquito diffusion only 
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(g) Percentage of infected cattle in the 
neighbors of Farm 3526 with 10000 in- 
fected eggs in Farm 3526 of Zavala 
county at the beginning considering 
mosquito diffusion only 




Tinne(Day) 

(h) Percentage of infected humans in 
the neighbors of Farm 3526 with 10000 
infected eggs in Farm 3526 of Zavala 
county at the beginning considering 
mosquito diffusion only 



Figure 10: The temporal dynamics of nodes in Uvalde area only with mosquito movement. The values are the aveage 
over 200 runs. 
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